Electronic correlations stabilize the antiferromagnetic Mott state in Cs 3 C6o 
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CS3C60 in the A15 structure is an antiferromagnet at ambient pressure in contrast with other 
superconducting trivalent fullerides. Superconductivity is recovered under pressure and reaches the 
highest critical temperature of the family. Comparing density-functional calculations with general- 
ized gradient approximation to the hybrid functional HSE, which includes a suitable component of 
exchange, we establish that the antiferromagnetic state of CS3C60 is not due to a Slater mechanism, 
and it is stabilized by electron correlation. HSE also reproduces the pressure- driven metalization. 
Our findings corroborate previous analyses suggesting that the properties of this compound can be 
understood as the result of the interplay between electron correlations and Jahn- Teller electron- 
phonon interaction. 

PACS numbers: 71.20.Tx,74.20.Pq,74.70.Wz 



Superconductivity in light-element materials entered 
in a new era after the recent observation of a critical 
temperature as high as 38K in CsaCeofH and the birth 
of a new family of organic superconductors based on 
aromatic molecules, as doped picene phenanthreneQ, 
coronenefH and dibenzopentacene[|3|. The unavoidable 
question about the "pairing glue" in the fullerides seems 
to have an answer at least for the standard members 
of the family, like K3C60 or RbsCgo- For these com- 
pounds there is convincing evidence of an electron- 
phonon driven superconductivity, in which Jahn- Teller- 
coupled intramolecular vibrations play the main role [6]. 
Nonetheless, the anomalous properties of "expanded" ful- 
lerides with large intermolecular distance, like ammo- 
niated compounds^ and most notably CS3C60 suggest 
that labeling these materials as standard superconduc- 
tors described by the Bardeen-Cooper-Schrieffer (BCS) 
theory is at least hazardous. 

The anomalies are particularly clear in CS3C60, in 
which the large Cs ions increase the lattice spacing in 
comparison with K- and Rb-doping. This compound 
is synthesized in two different crystal structure, the fee 
which is common to most doped fullerides and a A15 
structure, with a bipartite lattice^. A15 Cs3C6o at 
ambient pressure is an antiferromagnetic (AFM) insu- 
lator with Xjv = 46K and a spin- 1/2 moment at each 
molecule^, as opposed to K3C60 and Rb 3 C6o, which 
are superconductors in the same conditions. A super- 
conducting state is recovered under external pressure of 
around 4Kbar. Even more interestingly, the critical tem- 
perature has a bell-shaped behavior as a function of pres- 
sure, with a maximum of 38K at P = 7KBar. The prox- 
imity between superconductivity and antiferromagnetism 
and the existence of a maximum remind of the phase dia- 
gram of several exotic superconductors, ranging from the 
high-temperature copper-oxides to heavy fermion mate- 
rials, even if in the present case the bell-shaped curve 



is a function of lattice spacing instead of doping. How- 
ever, since the phase diagrams of these latter materials 
are believed to be dominated by strong electron-electron 
correlations and by the proximity to a Mott insulator, 
this similarity suggests that also the physics of CS3C60 
can be understood in the same framework. 

The role of electron-electron correlations has been 
underlined well before CS3C60 was known on the ba- 
sis of estimates of the Coulomb interaction and the 
bandwidth loL 11 1 In this light, the synthesis of CS3C60 
provides a system in which the correlation effects, already 
present in the standard fullerides, become predominant 
due to the reduced bandwidth associated to the larger 
distance between fullerene molecules. Once the role of 
correlations is realized, the AFM state of CS3C60 can 
be interpreted as a Mott insulator with subsequent or- 
dering of the Mott-localized electron spins. This pic- 
ture has indeed been drawn using Dynamical Mean-Field 
Theory (DMFT) [lj| for a three-band Hubbard model 
including the realistic electron-phonon interaction and 
the electron-electron correlation. These calculations have 
indeed predicted the phase diagram as a function of 
volume (or equivalently of the ratio U/W, where U is 
the screened Coulomb repulsion and W is the width of 
the conduction band), including the first-oder transition 
between a spin- 1/2 AFM and the superconductor (l3j . 
Within the same scenario, the correlation effects also 
influence the superconducting phase enhancing the or- 
der parameter, despite its s-wave nature 13l Il4j. The 
low-spin character of the insulator (the three electrons 
per molecule could also form a spin-3/2) is understood 
in terms of the Jahn- Teller interaction which favors the 



low-spin configuration 1J, |15 |. 

However, given the bipartite character of the A15 lat- 
tice, we can not rule out a Slater character of the AFM 
state. In the Slater picture the AFM ordering does not 
involve localized magnetic moment, but it arises as a 
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Fermi-surface instability of an underlying metallic state. 
A further confirmation of the correlated nature of the 
AFM state would also corroborate the idea that the same 
correlation effects arc crucial to understand the supercon- 
ducting state, while a Slater AFM would be more com- 
patible with a standard BCS mechanism for the super- 
conducting state. In order to further investigate the role 
of electron-electron interactions in the insulating state 
and to discriminate between a Mott and a Slater AFM 
state, in this work we tackle the physics of CS3C60 in the 
A15 structure from a complementary perspective, using 
density-functional theory (DFT) to investigate the ma- 
terial specific properties and its bandstructure. 

We anticipate that our main conclusion will be that 
local and semilocal functionals (Local density approxi- 
mation, LDA and generalized gradient approximation, 
GGA) are not able to properly describe the compound 
as an antiferromagnet, while using the hybrid functional 
introduced by Heyd, Scuseria and Ernzerhof (HSE)[16] 
an AFM state is obtained, as well as the pressure-driven 
metalization. These findings confirm that a pure Slater 
picture does not hold for the AFM state, and that the 
inclusion of correlation effects is a necessary ingredient to 
properly describe the properties of expanded fullerides. 

Before entering the details of our calculations, we 
briefly recall the main features of the electronic struc- 
ture. The bandstructure of doped fullerides reflects the 
molecular nature of the compound. The electronic bands 
correspond indeed to the molecular orbitals of an isolated 
Ceo molecule. The most relevant bands for the physics 
of doped fullerides are those originating from the three 
t\ u LUMOs of each Cqo molecule, which are expected 
to be populated by the three electrons donated by the 
alkali-metal atoms. The molecular nature of the solid re- 
flects into narrow bands, which lead to major influence 
of the sizeable intramolecular Coulomb repuls ion [id , Il7 1 
and also to a relevance of non-adiabatic effects [18| in the 
Jahn- Teller electron-phonon coupling which leads to the 
formation of Cooper pairs. 

The significant ratio U/W ~ 1.5 is a serious threat 
for density-functional calculations based on LDA and 
GGA approximations for the exchange-correlation po- 
tential. Therefore most studies devoted to the role of 
correlations in fullerides have resorted to many-body 
approaches to solve model Hamiltonians, ranging from 
Quantum Monte Carlo [l9j to Dynamical Mean-Field 

TheorytHQIIllill. 

On the other hand, the develop- 
ment of hybrid functionals including correlation effects is 
extending the realm of materials which can be handled 
within DFT providing us with a powerful tool which can 
be used to assess the relevance of correlations in a ma- 
terial. In this work we follow this route and we compare 
the results of GGA with the Heyd-Scuseria-Ernzerhof 
functional (HSE)0 as implemented in VASP 22!] in or- 
der to establish the role of electronic correlations in A15 

CS3C60- 
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FIG. 1. (Color online) Schematic view of CS3C60 unit cell 
with A15 structure showing the different orientation of the 
fulleride anions at the body center (B) and the origin of the 
coordinate system (A) and the position of the Cs cations 



The electronic structure of CS3C60 has been already 
investigated within LDA 17|, [23J . As a first step we com- 
plement the previous analyses using GGA in the Perdew- 
Burke-Ernzerhof (PBE)[2J| scheme. We use the projec- 



tor augmented wave scheme 25] and the cut-off for the 
plane- wave basis set was chosen as 400 eV and a 5x5x5 
mesh was used for the Brillouin-zone sampling. Calcula- 
tions are performed on the A15 structure of CS3C60 with 
Pm3n symmetry based on b.c.c. anion packing and re- 
fined at ambient pressure [l|. In this molecular packing 
the fulleride anion at the body center is rotated of 90 de- 
grees around the [100] direction with respect to the anion 
at the origin as schematically shown in Fig. [T] 

In Fig. [2] we report the density of states for the ti„ 
bands in the A15 structure obtained within PBE, which 
compare well with previous LDA calculations [l?], l23j |. 
The three bands arising from the ti u LUMO are well 
separated from the rest of the bandstructure. The band- 
width is around 0.55 eV, confirming the expectation of a 
narrow-band solid. The crystal field preserves the degen- 
eracy of the t\ u orbitals and three electrons per fullerene 
molecule are donated by the Cs atoms leading to a half- 
filled manifold which can, at least in principle, be unsta- 
ble toward antiferromagnetic (AFM) ordering with the 
two inequivalent sublattices associated to the two in- 
equivalent fullerene molecule. 

However, allowing for AFM ordering in the PBE 
scheme does not lead to the opening of an AFM gap. Im- 
posing an AFM magnetic structure with opposite spins 
in the two sublattices A and B (shown in Fig. [T}, only a 
metallic AFM solution can be stabilized. Moreover this 
state does not lower the energy with respect to the non- 
magnetic solution which remains the groundstate in this 
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FIG. 2. (Color online) Density of states obtained in PBE 
for A15 CS3C60 in a wide range of energies. The different 
bands arising from molecular orbitals are well visible and the 
tin bands are emphasized. The blue dashed line is the non- 
magnetic solution, while the red solid line is a metallic AFM 
solution 



approximation. Therefore, despite the bipartite nature 
of the A15 structure, nesting is not sufficient to stabilize 
the AFM insulator which is experimentally observed and 
a pure Slater scenario can be ruled out. 

In light of the significant value of the intramolecular 
Coulomb interaction, electron correlations can be invoked 
to explain the discrepancy between PBE calculations and 
experiments. In order to overcome this limitation, we re- 
sort to the HSE functional which introduces interaction 
effects through a suitable fraction of the exact exchange 
term. The nonmagnetic solution does not introduce sub- 
stantial differences with respect to PBE and only quan- 
titative effects modify the bandstructure (see the blue 
dashed line in Fig. [3J where we focus only on the t\ u 
bands) . The situation changes dramatically for the AFM 
solution, whose density of states is plotted as a solid line 
in Fig. [31 As opposed to PBE, an antiferromagnetic 
state is instead stabilized within HSE. The AFM state is 
an insulator with a finite gap A of around 0.25 eV and 
an energy gain with respect to the metal of around 0.3 
eV per unit cell. The ordered state has a magnetic mo- 
ment of ~ 1.4 [Ib per molecule, which is not far from 
the experimental value of ~ 1 [Lb M which corresponds 
to a spin 1/2 state. The discrepancy can be understood 
within the scenario of Refs. 13- lij], where the low-spin 
state has been attributed to the Jahn- Teller interaction 
which favors low-spin configurations in the Mott state. 
In our HSE calculations we did not relax the ionic posi- 
tion because this is computationally too demanding, and 
we expect this would lead to a Jahn- Teller energy gain 
and consequently to a reduction of the molecular spin. 

As the AFM state at ambient pressure is well repro- 
duced we also performed a series of calculations aiming 
at describing the experimental pressure-driven metalliza- 
tion. Of course our DFT calculations do not allow for 
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FIG. 3. (Color online) Density of states of the tin va- 
lence bands obtained with HSE functional in the nonmagnetic 
(dashed blue line) and antiferromagnetic (solid red line). The 
arrow pinpoints the AFM gap A 



superconductivity and describe the experimental super- 
conductor as a metal. The pressure-induced transition 
to bulk superconductivity takes place without change in 
crystal structure or symmetry [l|. The effect of pres- 
sure is therefore essentially a reduction of the lattice pa- 
rameter, which can be easily taken into account within 
DFT. We performed calculations for the lattice param- 
eters corresponding to 0, 5 and 10 kbar as determined 
in Ref. [l[ without relaxing ionic positions. The re- 
duced lattice spacing leads to a larger overlap between 
the ti M orbitals of different molecules and consequently 
to a larger bandwidth, which becomes 0.645 at 5 kbar 
and 0.685 at 10 kbar. The increased bandwidth leads 
to a reduced effect of correlations, which in turn modi- 
fies the energetic balance between the nonmagnetic and 
AFM states. At 5kbar the two solutions have the same 
energy within our accuracy, while at 10 kbar the metal is 
stable with an energetic gain of more than 0.3 eV per unit 
cell. This evolution remarkably mirrors the experimen- 
tal findings and it strikingly confirms how a significant 
degree of correlation is necessary to stabilize the AFM in- 
sulator. Interestingly, at 5 kbar the magnetic moment of 
the magnetic solution is still 1.34 /is, only slightly below 
the ambient pressure value, despite the system is at the 
boundary between the two phases. This points towards 
a strongly first-order transition between an AFM with a 
finite magnetic moment and a metal which is expected 
to become superconducting. 

HSE is thus able to reproduce both the insulating be- 
havior at ambient pressure and the metallic phase at 
finite pressure. Moreover, it provides an estimate of 
the transition which is close to the experimental result. 
Again, the discrepancy between our calculations and ex- 
periments can be attributed to the lack of ionic relaxation 
and the consequent underestimate of the Jahn- Teller en- 
ergy gain. Including the Jahn- Teller interaction is in- 
deed expected to reduce the magnetic moment and to 
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slightly increases the energy of the AFM state. As a con- 
sequence, a calculation including the full Jahn- Teller dis- 
tortion would provide a smaller critical pressure for the 
metalization. The ability of HSE to essentially reproduce 
the evolution of Cs3Cgo under pressure is a strong evi- 
dence of the accuracy of this approach to address strongly 
correlated molecular solids. This opens the path for a 
number of extensions which can bridge the gap between 
DFT calculations and, for example, DMFT solutions of 
a three-band Hubbard model. 

The success of HSE for CS3C60 can be used as an indi- 
rect indication of the accuracy of the same approach also 
for picene and for other organic superconductors, which 
all share similar values of ?7/W[17|. A HSE calculation 
for potassium-doped picene indeed confirms the relevance 
of correlations [26j. HSE has also been used to propose a 
positive role of electron correlations in determining the 
critical temperature of Ba^ K x B1O3 and -HfNCl 27]. 

In this paper we have demonstrated that HSE hybrid 
functional is able to reproduce the evolution of the prop- 
erties of A15 CS3C60 as a function of pressure, while PBE 
is not. In particular, HSE finds an AFM state at ambi- 
ent pressure and an insulator-to-metal transition under 
external pressure, while PBE always gives a metallic so- 
lution. The calculations give a magnetic moment which 
is slightly larger than the experimental value, a difference 
which can be ascribed to the lack of ionic relaxation in 
our calculations. Only the full distortion of the molecule 
would allow for a full Jahn- Teller energy gain and a spin- 
1/2 local moment. Our results strengthen the idea that 
electronic correlations are at the basis of the properties 
of Cs3C6o- This provides an indirect confirmation of 
the strongly correlated superconductivy scenario [HQ. 
The similarity between the present results and previous 
calculations for potassium-doped picene suggest that the 
aromatic superconductors may share many similarities 
with doped fullerides. A further experimental test of 
the Mott nature of the AFM state can be provided by 
spectroscopic methods. It has indeed been shown theo- 
retically that a three-dimensional Mott antifcrromagnet 
should display a peculiar series of spin-polaron satellites 
both in the photoemission spectra[28J and in the optical 
conductivity, in which equally spaced spin-polaron peaks 
are predicted to appear above the optical gap 29] . 
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